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^ ■ Abstract 

ps| I The multiplicity distributions at high energy e~^e annihilation are described well 

within the Two Stage Model in the region from 14 to 189 GeV. Energy dependence 
^ ■ of parameters of this model gives the dynamic picture of the parton stage and the 

, stage of hadronization. It is shown that oscillations in sign of the ratio of factorial 

I cumulant moments over factorial moments of increasing order can be confirmed by 

' this model. 

o 

O , 

^ ■ 1 Introduction 

Multiparticle production (MP) is one of the most important topics in high energy physics. 
Using MP we can get more information about the nature of strong interactions and un- 
derstand deeper the structure of matter. Over the last few years many thorough reviews 
devoted to MP have been done (Ij. Modern accelerators have made it possible to study 
rS I MP more extensively and in detail. Developing theory of high energy physics quantum 

I chromodynamics (QCD) |2] and a lot of phenomenologic models are tested by the process 
of MP. 

Multiparticle processes begin at high energy. Among all producing particles we can 
observe a lot of hadrons. On one hand we want to know high energy physics, but on 
the other hand the increase of the inelastic channels makes it difficult to describe this 
process with customary methods. The situation concerning the history of thermodynamics 
developing and statistical physics is much the same. Analysis of MP process is carried out 
using of statistical methods because the number of secondaries in e+e~ annihilation is 
large (more than 60) [3 . The consideration of MP begins from the behavior of charged 
multiplicity. 
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As it is generally known the multiplicity is the number of secondaries n in process of 
MP: A + B ai + a2 + . . . + The multiplicity distribution (MD) Pn is the ratio 
of cross-sections cr„ to a = '^CTn'- Pn = cTn/o"- This quantity has the following meaning: 

n 

the probability of producing of n charged particles in this process. We can also construct 
quantities such as mean values, moments of MD, can study correlations and so on. 

Investigation of MP has led to discovery of jets. Jets phenomena can be studied in 
all processes, where energetic partons are produced. The most common ones are in e~^e~ 
annihilation, deep inelastic scattering of e, fi or u on nucleons and hadron- hadron scatter- 
ing, involving high-px particles in final state. Let us consider e~^e~ -annihilation at high 
energy. This process is one of the most suitable for the study of MP. In accordance with 
QCD it can be realized through the production of 7 or Z°-boson into two quarks: 

e+e- (Z77) ^ gg (1) 

Perturbative QCD can describe the process of fission partons (quarks and gluons) at 
high energy, because the strong coupling as is small at that energy. This stage can be 
called as the stage of cascade. After hard fission, when partons have not high energy, 
they must be changed into hadrons, which we can observe. On this stage we shouldn't 
apply perturbative QCD. Therefore phenomenological models are used for description of 
hadronization (transformation of quarks and gluons into hadrons) in this case. 

The description of the stage quark-gluon cascade by means of perturbative QCD was 
applied in [H E] . Certain features of the predictions at the parton level are expected to be 
insensitive to details of the hadronization mechanism. They were tested directly by using 
hadron distributions jH^. 

The e"'"e~-reaction is simple for analysis, as the produced state is pure qq. It is usually 
difficult to determine the quark species on event-by-event basis. The experimental results 
are averaged over the quark type. Because of confinement the produced quark and gluons 
fragment into jets of observable hadrons. 

The hadronization models are more phenomenological and are built on the experience 
gained from the study of low-px hadron collisions.lt is usually considered that the pro- 
ducing of hadrons from partons is universal process. 

2 Two Stage Model 

Parton spectra in QCD quark and gluon jets were studied by Konishi K., Ukawa A. and 
Veneciano G.W. Working at the leading logarithm approximation and avoiding IR diver- 
gences by considering finite x, the probabilistic nature of the problem has been established 

la 

At the studying of MP at high energy we used idea of A. Giovannini [7J for description 
of quark-gluon jets as Markov branching processes. Giovannini proposed to interpret the 
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natural QCD evolution parameter Y 



^ ln[l + a,61n(^)], 



2nb 



(2) 



where 27ib = ^{llNc — 2Nf) for a theory with Nc colours and Nf flavours, as the thickness 
of the jets and their development as Markov process. 
Three elementary processes contribute into QCD jets: 

(1) gluon fission; 

(2) quark bremsstrahlung; 

(3) quark pair creation. 

Let AAy be the probability that gluon in the infinitesimal interval AY will convert 
into two gluons, AAF be the probability that quark will radiate a gluon, and BAY be the 
probability that a quark-antiquark pair will be created from a gluon. A, A, B are assumed 
to be Y-independent constants and each individual parton acts independently from others, 
always with the same infinitesimal probability. 

Let us define the probability that parton will be transformed into m gluons over a jet 
of Y in thickness and call it P^{Y). The probability generating function for a parton jet 
will be 



Q^{z-Y) = Y,P^{Yy. 



(3) 



m=0 



A.Giovannini constructed system of differential equations and obtained explicit solu- 
tions of MD for a parton jet in particular case 5 = (process of quark pair creation is 
absent) In the common case -B 7^ MD are similar to particular one (Tj. 

For quark jet explicit solutions are given [7j 



-AY 



PmiY) 



/i(/i+l)...(/^ + m-l) ^_^y _ ^-AY^m 



ml 



(4) 



where /x = ^ . Futher the average gluon multiplicity is m = /i(e'^^ — 1) and the normalized 
exclusive cross section for producing m gluons from quark is 



0-^ 



O'tot 

+ 1) . . . (/i + m - 1) 
ml 



Prn{Y) 



m 



m + /i 



m + jj 



The generating function ^ will be given by 

00 

Q'^(^,Y) = ^z™P„(Y) = 



m=0 



-AY 



l-z(l-e- 



-AY\ 



(5) 
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Eq.(jH) is Polya-Egenberger distribution, where is non-integer. 

In Two Stage Model [S] we took (jH) for description of cascade stage and added su- 
pernarrow binomial distribution for hadronization stage. We chose it based ourselves on 
analysis of experimental data in e+e~- annihilation lower 9 GeV. Second correlation mo- 
ments were negative at this energy. The choise of such distributions was the only could 
describe experiment. 

We suppose that hypothesis of soft colourless is right. We add stage of hadronization 
to parton stage with aid of it's factorization. MD in this process can be written 



(7) 



where is MD for partons (0]), P^{m,s) - MD for hadrons produced from m partons 
on the stage of hadronization. Futher we will use instead of parameter Y CM(center of 
masses) energy ^/s. 

In accordance with TSM the stage of hard fission of partons is described by negative 
binomial distribution (NBD) for quark jet 



kp{kp + 1) . . .{kp + m - 1) 



m 



m + kr 



kn 



k„ + m 



where kp = A/ A, m = '^rriP^. We neglect process (3) quark pair production {B = 0). 

m 

Two quarks fracture to partons independently of one another. Total MD of two quarks 
is equal to ((Zj) too. Parameters kp and m of MD for two joint quark-antiguark jets are 
doubled, but we use that designations. 

and generating function for MD Q^{s, z) are 



1 9"" 
m! dz^ 



QlX^.z) 



X _l_ (^X 
kp 



MD of hadrons formed from parton are described in form |Sj 



(9) 
(10) 



C 



n:. 



N„ 



{C^ - binomial coefficient) with generating function 



1 + 



n' 



(12) 
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where and Np {p = q, g) have meaning of average multiphcity and maximum secondaries 
of hadrons are formed from parton on the stage of hadronization. MD of hadrons in e'^e~ 
annihilation are determined by convolution of two stages 

°° an 

Pn{s) = Y.P^^{Q^?'--l=.. (13) 
m=o 

where 2 + m is total number of partons (two quarks and m gluons) . 

Further we do the following simplification for the second stage: pa considering 
that probabilities of formation of hadron from quark or gluon are equal. We introduce 
parameter a = jf- for distinguishing between hadron jets, created from quark or gluon on 

the second stage. We also make simplification for designation = Ng, = n^. Then we 
get 

<3f=(i + |(^-i))", 

Introducing in p3|) expressions (jH)), (fT^ and differentiating on z we obtain MD of hadrons 
in the process of e^e~ annihilation in TSM 

^n(3) = $^P^C2\_)^ - ■ (14) 

m=0 V / V / 

For comparing with experimental data the normalized factor Q was introduced into (fT!?|l 
and a number of gluons in the sum was restricted by Mg - maximal number of possible 
gluons created on the first stage 

m=0 V / V / 

The results of comparison of model expression (fT3j) with experimental data J2] are 
represented in Table 1 (parameters of two stages) and on Figures (1)-(18). We can see that 
MD in TSM (solid curve) are describing well the experimental e"'"e~-data (black square ■) 
from 14 to 189 GeV. Summing up is limited to Mg equal 20 — 22 for energies to 61.4GeV 
and 33 — 41 above, because further increase does not change x^- 

The description of experimental MD ^2] by (fTH|l gives ~ 1 for the most energies 
with the expection of 34.8GeV and 55GeV (more high x^), Ni:,p = 14. Calculated MD are 
giving small deviations from experimental data at 34.8GeV, but for almost all multiplicities 
(Fig. 3). Probably it is connected with small statistics of events. 

At 55GeV they give big deviations, but only to the right side of central region (Fig. 7). 
It can be described as energy suppression of formation of hard gluons connected with 
appearance of heavy qq pair or heavy hadrons. Energy behavior of total cross section of 
e'^e^ hadrons P", is changed to a sharp rise from here. It also may be connected with 
the reason of narrowing of distribution. 
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Table 1. Parameters of TSM. 
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3 Dynamics of multiparticle production 

We will analyse dynamics of MP which corresponds to values of parameters of TSM 
(Table 1). We will begin from a cascade stage. This stage is described by two parameters: 
m and kp. 

The average multiplicity of gluons m formed on fission stage has tendency to rise. It 
is changed from ~ 0.1 at 14GeV to ~ 20 at 183 Gev. But we can see certain insignificant 
deviation from this direction at ^/s=50 — 61.4GeV, at 183 GeV. It foUowes from QCD 
that in particular case [B = 0) the parameter kp that equal to ratio 2 A/ A — > 1. Values kp 
are changed insignificantly. They are remained ~ 10 at almost all energies. There is some 
physical senses of this parameter. One of the most interesting from them is temperature T 
|10j : T = k~^. From thermodynamical models we can also obtain the following connection 

k;' = To + l/cE, (16) 

where Tq is the temperature of system before interaction, c - thermal capacity, E - energy 
spended on creationg new particles JU]- In this sense we can make assumptions: temper- 
ature of parton system with developed cascad are lowest at 14 and 55GeV than at the 
others. The highest temperatures are reached on first stage at 50, 52, 61.4 and 183Ge\^. 

The interesting picture of hadronization is discovered in conformity with parameters of 
second stage Nq, and a . The parameter Nq determines maximum number of hadrons, 
which can be formed from quark on this stage. In TSM (Table 1) it takes different values 
from 4 to 55. We can't reveal steady energy rise or fall for it. Big A^^ point to predominance 
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of hadrons formed from quark jets, small A'^, point to essential contribution gluon jets in 
the hadron multiplicity. More probably that this parameter remaines constant and ~ 16 
with small deviations. 

Next parameter has meaning of mean hadron multiplicity from quark on second 
stage. We can see the tendency to weak rise with big scatter. Such behavior of parameter 
may be connected with the growth of spectrum of mass hadron states (appearance of new 
mass states with increase of energy). The average value of is about 5 — 6 in the research 
region. 

Parameter a was introduced for comparision quark and gluon jets. It is almost constant 
and equal to 0.2 with some deviations. If we know a then we can determine analogous 
parameters A^^^ = aNg and = an^ for gluon jet. It is interesting that these parameters 
remain constant without considerable deviations: Ng ^ 3 and ~ 1 (Figures 19 and 20). 
From this result we can affirm about universality of hadronization. 

The fact that a < 1 says that hadronization of gluon jets are more soft than quark one. 
The simplest explanation to this phenomenon is the fact that a quark takes away more 
conciderable energy than gluon. 

The ratio determines the probability of formation of hadron from parton. It is 
increased from ~ .1 to ~.5 with rise of energy from 14 to 43.6GeV, then we have big 
variations in region 50 — 61.4 GeV (.23 — .69) and the ratio is almost constant at higher 
energies (~ .4) (Figure 21). It should be noted that the small probabilities are realized at 
55GeV (in the region 50 - 61.4GeV) and at 183GeV. 

The normalized factor Q remains constant and is equal to 2. 

4 Oscillation of moments in MD 

It was shown recently ^H] that the ratio of factorial cumulative moments over factorial mo- 
ments changes sign as a function of order. We can use MD formed in TSM for explanation 
of this phenomenon. 

The factorial moments can be obtained from MD P„ through the relation 

oo 

Fg = J2n{n-l)...{n-q + 1)P„, (17) 

n=q 

and factorial cumulative moments are found from expression 

q-l 

K, = F,-Y,CUKq^^Pi- (18) 

i=l 

The ratio of their quantities is 

H, = KJF,. (19) 
We can use the generating function for MD of hadrons (|Tl|l in e^e~ annihilation G(z) 

m=0 
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= Q^{Q';{z)rQl{z). 
We are calculating Fq and Kq in TSM using fl^ 

1 <9«G 
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z=Q 



' W{s) dzi 
The expression (120]) for G(z) after taking a logarithm 
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lnG'(s,^) = -fcpln[l + 7^(1 - gf )] + 21ngf 
and the expansion to series in power on will be 
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and using (j22j) we obtain 
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where n(s) is the average multiplicity hadrons in process ([T}. It is possible to find Fg 

m=0 ^ ^ 



with Pm equal (jHl). 

The sought-for expression for Hg will be 
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where Qi is the normahzed factor. The comparison with experimental data jlH] shows that 
()26|) describes the ratio of factorial moments (Figures 22-39). It is seen minimum at q=5. 

In the region before Hq may oscillate in sign only with the period equal 2, changed 
sign with parity q. At more high energies the period is increased to 4 and higher. It can 
be explained by influence of hadronization. Values Kg (as well Hg) may change sign only 
owing to second summand in (j^^ . More devo loped cascad of partons with hadronization 
comes to big period of oscillations of sign. 

The immediate calculations Hg based on (fT7jl - (fT^ with using MD (P3j) gives very good 
description of the oscillation value of Hg {x^ ^ 2). Significant oscillations start near region 
producing of Z° and can be explained by non-integer values of parameters of hadronization 
Ng and Ng = aNg or by convolution of wide (for parton jets) and narrow (for hadron jets 
on second stage) MD. 

5 Conclusions 

It is shown that TSM does not contradict to the experimental data on MD and the oscil- 
lations ratio of factorial moments. TSM offeres concreted physical picture of multiplicity 
production in high energy e'^e~ annihilation. 
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Figure 1: MD at 14GeV. Figure 2: MD at 22GeV. Figure 3: MD at 34.8GeV 
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Figure 4: MD at 43.6GeV. Figure 5: MD at 50GeV. Figure 6: MD at 52GeV. 
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Figure 7: MD at 55GeV. Figure 8: MD at 56GeV. Figure 9: MD at 57GeV. 
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Figure 17: Parameter A'^ = aNq. 



Figure 18: Parameter = amJl 
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Figure 19: Ratio 
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Figure 20: Hg at 14GeV. Figure 21: Hq at 22GeV. Figure 22: Hg at 34.8GeV 
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Figure 23: Hg at 43.6GeV. Figure 24: Hg at 50GeV. Figure 25: Hg at 52GeV. 
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Figure 26: Hg at 55GeV. Figure 27: Hg at 56GeV. Figure 28: Hg at 57GeV. 
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Figure 29: Hq at 60GeV. Figure 30: Hq at 60.8GeV. Figure 31: Hq at 61.4GeV. 
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Figure 32: Hq at 91.4GeV. 
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Figure 33: Hq at 172GeV. Figure 34: Hg at 183GeV. Figure 35: Hq at 189GeV. 
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